LMU-ASC 09/07 



The high momentum behavior of a quark wake 
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Using the AdS/CFT duality, we analytically evaluate the high momentum (or 
short distance) behavior of the color field strength due to a moving quark in an 
N = 4 SYM plasma. We find a fireball-like behavior in the near quark vicinity as 
predicted in the literature on general grounds, in the context of heavy ion collisions. 
Our approach to analytically solving the problem is based on a WKB approximation 
which may be extended to other setups as well. 



I. INTRODUCTION, SUMMARY AND DISCUSSION 



X 



Recent results from the Relativistic Heavy Ion Collider (RHIC) at the Brookhaven Na- 
tional Laboratory [l[ 4^ (see j|J for a review) have lead to an increasing interest in 
the details of the strong interactions at temperatures not too far above the deconfmement 
transition Tqcd ~ 170 MeV. At these temperatures the gauge coupling g\ M is, apparently, 
large. One approach towards understanding this type of strongly interacting phenomenon 
is through the AdS / CFT duality, where one considers the high temperature phase of large 
N SU(N) M = 4 Super Yang Mills (SYM) theory as a first step towards understanding 
Quantum Chromodynamics (QCD). While it is not completely clear if, how or when these 
two theories may be likened 6], an ad hoc comparison of some of their features gives re- 
markably similar results. This includes an analysis of the energy density at strong coupling, 



the viscosity to entropy density ratio 
quarks and mesons [ill . Il2l 
the jet quenching parameter 




10 L the diffusion or friction coefficient for 



an analysis of meson melting 



29j , the screening length 
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2l|, 



and 
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directional emission [32|, |33j, |34| . See also |35|, |36|, [37|, [38|, [39( for other properties of M = 4 
SYM which were discussed in the literature in this context. 

Being motivated by the relative success of the above analyses, we have considered the 
directional emission of colored fields due to the motion of an infinitely heavy quark in an 
M = 4 SYM plasma. As the quark moves through the plasma it emits gluons which may be 
characterized by their directionality and energy. A numerical analysis of such a gluon wake 
was initiated in {32], and continued in 33] and 34] where a comparison to RHIC results 
was elaborated on. In this note we continue this investigation and, using a WKB 1 based 
approach, obtain an analytic result for the directional emission of colored fields at large 
momentum scales. 

The magnitude of the field strength in response to the motion of the quark, TrF 2 , may be 
studied by considering the related, Fourier transformed (hatted) quantity (0$(A;_, k±)) = 
— — Tr (F 2 + derivative ) to be defined in detail in section [III Here fc_ is the momentum 

$Y M \ terms / 

component conjugate to x_ = x — vt, where x is the direction of motion of the quark and 
v is its velocity, and k± is the momentum conjugate to the directions perpendicular to the 
motion of the quark. We find 



YM 



1 



+ i 
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\/ AyM (1 



(1 - V 2 )k 2 _ + k\ 

~) nvk_T 2 



0(T 4 /((l-v 2 )k 2 _ + k 2 ± ) 2 ). (1) 



3 (1 - v 2 )k 2 _ + k\ 

In equation ([!]) T denotes the temperature of the plasma and Xym the 't Hooft coupling 
constant of the SYM theory. The first term is identical to the near field contribution of a 
moving quark in a conformal theory at zero temperature {32], 41], while the second term is 
a purely dissipative effect. We shall concentrate on the latter. 

A contour plot of the dissipative component of (0^(k_,k_ s _)} is given in figured! The 
directional behavior of this quantity may be analyzed by considering the gradient of (0$) 

d 

(where the 'd' implies that we are considering the sub-leading dissipative contribution). In 



the k± direction the gradient of (0$) 
extremized at 

tan# 



is extremized at k 1 



k 



0. In the k- direction it is 



(2) 



1 Wentzel-Kramers-Brillouin- Jeffreys, see for example 40 1 
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velocity 


0.75 


0.90 


0.95 


0.99 


Numerical estimation of 9 [32] (k-Zo = 40) 


0.58 


0.41 


0.30 


0.17 


Analytic approximation of 9 (large k^z^) 


0.58 


0.41 


0.30 


0.14 



TABLE I: We compare our results for the angle of the high momentum ridges given in equation 
([2]) to the numerical results presented in [32:]. As the velocity approaches the speed of light, the 
regime of validity of our approximation shifts towards larger values of the momenta (see section [TT] 
for details). The last column shows that for v = 0.99, the quoted value of 9 ~ 0.17 at = 40 32J 



will decrease as ^ increases beyond 40, until it will stabilize at an angle of 9 = 0.14. 



This directional behavior may be graphically enhanced by multiplying the momentum space 
density (0$) by k±; since k±(0<s>) is of order k° ~ 1 it will not decay at large values of 

d d 

the momenta. The gradient in the k- direction is unaffected by this amplification, though 
the gradient in the k± direction is modified reaching a minimum at k± = 0. Topographically, 
one finds a 'ridge' in momentum space at an angle 9 given by equation ([2]). A comparison 
of the angle 9 to the numerical value found in [32J is given in table [B 

The momentum space distribution ([I]) may be Fourier transformed to real space to obtain 
the near quark value of the magnitude of the color field strength. One finds that 

y/X YM (l - V 2 )vX„T 2 



(0.) 



YM\ 



,2\2 



+ 0(T 4 (r 2 (l-v 2 )+x 2 _) 2 ) . (3) 



16tt 2 (r 2 (l - v 2 ) + x 2 _f 12 (r 2 (l - v 2 ) + x 2 _f' 2 

Considering again the sub-leading dissipative term, we find that the gradient in the r direc- 
tion is maximized at r = 0, while in the X- direction it is maximized at an angle 

V2 



tancu = — 



(4) 



see figure [2j 

As discussed in 



32 



34| , it is interesting to consider this result in light of the suppression of 
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451 ]: one 



the back-to-back di-hadron correlations found in the RHIC experiment [42|, 
expects that jets created as a result of parton collision will be emitted from the quark gluon 
plasma (QGP) at an angle of 180° to each other. Instead, one finds that when considering 



jets with transverse momenta in the range 2.5GeV/c < pt < 4.0GeV/c, one of the jets is 



missing or seems to be deflected by an angle of around 1 radian 



46 



471 ] . This behavior is 
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FIG. 1: Contour plots of the momentum space distribution of the leading dissipative behavior of 

the dimensionless quantity ~ T v / I^7 2glr M ^ + • ■ ■) ^ n res P onse t° a 

quark moving at a constant velocity of v = 0.8. The horizontal axis corresponds to the momentum 
conjugate to the comoving coordinate of the quark, while the vertical one corresponds to 
momentum conjugate to the direction transverse to the quark motion, k±. The momentum is 
measured in units of temperature (denoted [T] in the plots). Positive imaginary values are shaded 
dark while imaginary values closer to zero correspond to lighter shading. The origin of momentum 
space has been colored signifying that this region can not be accessed by our high momentum 
approximation. Figure (a) depicts the leading dissipative contribution to (0$). The red line 
follows the extremum of the gradient of (0$) in the k- direction. This directional behavior may 
be graphically enhanced (figure (b)) by multiplying (0$) by k±. 
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<0^ (x_ , r) > 



40 



x_ [1/T] 



FIG. 2: Contour plots of the leading dissipative behavior of the dimensionless quantity 

■=a A (0<b(x — vt.r)) = \ — =7rk — (TrF 2 + ...") in response to a quark moving at a con- 

stant velocity of v = 0.8. The horizontal axis corresponds to the comoving coordinate of the quark 
x — vt, while the vertical one corresponds to the direction transverse to the quark motion, r. Dis- 
tances are measured in units of inverse temperature (denoted [1/T] in the plot). Large negative 
values are shaded in dark while values closer to zero correspond to lighter shading. The red line 
marks the extremum of the gradient in the x_ direction. 



interpreted as an indication that the original collision occurred near the surface of the QGP, 
so that while one of the partons managed to exit the QGP, the other had to traverse it and 
as a result has been deflected or stopped. 

One attempt at explaining this deflection angle is related to the sonic boom created when 
the quark traverses the QGP at a speed higher than the speed of sound in the plasma 48]. 
In this case the sound wave created far from the moving parton will be directed at a Mach 
angle relative to the direction of motion of the parton which turns out to be approximately 
one radian. Indeed, in [34J] it was shown that for a quark moving in an J\f = 4 SYM plasma 
the low momentum contribution to the energy density has a preferred directional behavior 
at an angle of cos 2 9 ~ the Mach angle for a quark in a conformal QGP. 

The analysis of the sound wave resulting from the motion of the parton 48] takes place in 
a region far from the jet where perturbations and gradients are small and the hydrodynamic 
approximation is valid. In contrast, in the near parton region the estimated energy deposition 
of the parton to the medium is large, probably due to strong dissipative effects. This implies 
that the jet is surrounded by a "fireball" where the hydrodynamic approximation is not valid. 
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It is interesting to relate this fireball to the behavior of the near field of the parton given 
in equation (j3J) and depicted in figure [2j The directionality of (G$) (equation (jlj)) indicates 
that this fireball exhibits Lorentz-like contraction at relativistic velocities. 

In order to proceed with an analysis of the short distance behavior of the quark wake, one 
needs the thermal expectation value of the energy momentum tensor which was obtained 
numerically in [34J]. The energy momentum tensor will give us a better indication of how 
energy is distributed by the quark as it travels through the QGP than the field strength. This 
problem is somewhat more complicated due to the couplings between various components 
of the (dual) metric tensor's equations of motion. We plan to pursue this in the future. 

Finally, we point out that one should be somewhat cautious when relating these high 
momentum results, which are relevant to Af = 4 SYM, to real world QCD, especially when 
considering the high momentum ridges which correspond to a weakly coupled regime of the 
theory. 



II. EVALUATION OF (0$) 



In this section we evaluate the thermal expectation value of the operator Oq> 



Tr F 2 + den 



vativc 
terms 



in the presence of a moving quark in a thermal plasma. 



Using the AdS / CFT dictionary 



49 
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51 



521 ]. the motion of a quark moving at constant 



velocity in the x direction is given by a string in an AdS-Schwarzschild (AdS-SS) geometry 
whose endpoint is moving with constant velocity v along the asymptotically AdS boundary. 
The profile of such a string can be found by minimizing the Nambu-Goto action, 



S 



NG 



/e 2 v/j — h\drda. 



(5) 



In equation flS} we used h a p = G )lu d a X fJ 'dpX 1 ' with G the AdS-SS metric tensor which may 
be read off from the line element 

I - n ( 7\r\i 2 4- dr 2 4- ' 

y 2 



ds' 



L 2 dz 2 

-g(z)dt 2 + dx 2 + dx 2 + — 
:- \ ' ft .'/("-) 



(6) 



with g 



~0 



In this notation, the radial coordinate z lies in the range < z < Zq, 



with zq the location of the 
Solving for A M 



Dlack hole horizon. 
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121 ]. one finds that the string configuration is given by X 
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(t,£(z) + vt,0,0,z) with 



= ^ (in + 2 arctan - ) . ( 7 ) 

4 V + t ^ 



(See also llq for some more elaborate configurations which include the presence of other 
fields or a chemical potential.) 

The AdS/CFT duality then tells us that the thermal expectation value of the CFT 
operator (0$) is given by 

<°»> = -isk53>*- (8) 

where $ is the dilaton field in AdSs which is sourced by the moving string and G§ is the five 
dimensional gravitational constant (see equation (|T0|) below). In equation (jSJ) one should 
use 



N 



2 



, >n , ~™ = z o and v^m = — , (9) 
ZUr5 Til a 

with T the temperature, N the number of colors and Xym the large N 't Hooft coupling 
constant. 

The equation of motion for the dilaton ,$(x — vt, x%, X3), may be obtained by minimizing 
the action 

S = Sdilaton + SjVG 

with 

Sdilaton = -T^TT I V^G\kd$) 2 d 5 X. (10) 

lb7rG 5 J 2 

The resulting equation of motion for $ = 4G °!yjz^T / ^e~ tkp d?k with k = (A;_ , fc 2 , ^3) and 
p = (x - vt,x 2 ,x 3 ) is 

(>d 2 {z-*g(z)d M ) - (V - ^ 2 ) ) $ = *f*-««. (11) 

In general one finds that when expanding in small z, $ = M + iVz 2 — \Pz 3 + Qz 4 + 0(z 5 ). 
The constant M is related by the AdS / CFT duality to a deformation of the Lagrangian (a 
source term) and in the setup we are considering should be set to zero. This implies also 
N = 0. Using ([9]) and assuming that P depends analytically 2 on the momenta k% and 
fc 3 , equation (jSJ) becomes 

(Ch) = --y/\ YM (l-v*)Q. (12) 



7T 



2 Divergences which are analytic in the momenta a re, after a Fourier transform, contact terms and may be 
ignored (see 5l|, [53|, l54[ for general remarks and 55|, [56[ for some recent applications) . 
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The solution to equation (TTTT) has been evaluated numerically in 32] and analytically for 



small momenta in 



34|. In what follows we shall solve it for k > T. 



Before doing so, we give a physical explanation of our result. This will be followed by 
a more elaborate proof. We start by noting that as we take the dimensionless quantity 



a/ k 2 — v 2 k?LZo to be very large (or equivalently g — ► 1 or T — > 0), we obtain the equation 
of motion for a string hanging straight down in empty AdS space. The solution to this 
equation is uniquely specified by requiring that it be normalizable and source free (i.e., that 
it vanishes at the asymptotically AdS boundary). One may guess 

^ = W 2 fa 1 *^ - L ^ z)) + 2L i(^)) > ( 13 ) 

where we have introduced k 2 = k 2 — v 2 k 2 _, and Li is a modified Struve function of order i 
[57I . From this one can read off P = 1 and Q = j^irk so that using (TT21) one has 

(Oz) = -^\ YM (l-v 2 )±-k. 

lo 

Next, we note that expanding the equations of motion in inverse powers of 
\Jk 2 — v 2 k 2 _z = kz , the leading contributions to the left hand side of ffTTj) are of order 
(kz )~ 4 . Since ze~ lk -^ = z + % -^ r z A + O(kz )~ A , the leading contribution to the right hand 
side of (11 II) is of order (kzo)~ 2 . At this order the only solution which does not diverge 
exponentially in the interior and is not sourced is given by the combination 

$ = ™ (kz{I 2 {kz) - L (kz)) + 2Li(te)) - ^- 2 z\ 
lk A v J Sk^ZQ 

which implies that the sub-leading terms in (0$) are precisely those given in equation (CQ). 

Physically we have the following picture: momenta which are large relative to the tem- 
perature are associated with a region of spacetime which is close to the asymptotically AdS 
boundary. Therefore, if one probes large enough momenta, he or she will not notice that the 
geometry deviates from empty AdS and an analysis of the moving string will be identical to 
the one in AdS space. Finite temperature corrections are then encoded in the shape of the 
stringy source near the AdS boundary and not in the spacetime geometry. 

In the above discussion there is one issue that we have glanced over and that should be 
treated with more care. We have expanded the equation of motion in a power series in (/czo) 1 
although the normalizability condition on $ is imposed at the horizon, located at z = zq. 
This proves to be a non perturbative statement; once we expand the function $(^0, k; z) in 



9 



a power series in (fc^o)" 1 , we can no longer correctly evaluate &(zq, k; zq) without resuming 
the power series. Nonetheless, despite the somewhat hand-waving nature of this argument 
we shall see in what follows that the final result is correct. 

To solve (JTTJ) in the large kz limit, we find it useful to switch to the dimensionless 
variables 3 Z = kz, Zq = kzo, and define a = v^-. Redefining k 2 & = (f> one obtains the 
equation of motion 



Z 3 d z (Z- 3 g(Z)d z ) -1 + a 2 



g{z) 



- 1 



— e 
k 



-ik-£(Z/k) 



(14) 



Our approach to solve this equation is to first solve the appropriate homogeneous equa- 
tion, once with boundary conditions such that the solution has no source, 0_, and once 
with boundary conditions such that the solution has no outgoing modes at the horizon 4 , 
<f> + . From these solutions, the Greens function, G(Z,Z'), for equation ([141 or (TTT1) may be 
constructed, 



G(Z, Z>) 




where W(<f>-, <j>+) = <p' + 4>- - <j>'-' 
inhomogeneous equation for $ is then given by 



-(Z')MZ) z>z> 
(z')^(z) z<z>, 

is the Wronskian of </>+ and 6_. 



(15) 

The solution to the 



$(z) = k' 



kz 



G(z~k,Z')Z'e- ik ^ {z ' rk) dZ'. 



This method has been used in (4l| to obtain the spacetime form of f|T3|) . (0$) — 16pi . 
a motionless string (recall that p = (x-,X2,Xs)). 



for 



To find the homogeneous solutions to (IT41) we shall use the WKB approximation. This 



method has proved to work well when computing t 
to 59], or for more recent work [60] to 



41 



61 



re mass spectrum of glueballs (compare 



62|). 



In a quantum mechanical context the WKB approximation is valid provided one is con- 
sidering a high enough momentum or, conversely, a flat enough potential relative to the 



For comparison, the authors of 



32 



used y = j-, K\ — k_Zo and K± = k±z , so that Kf(l 



v 



D 



K 



K\ = Zq. The numerical solution to Bk of 32| or Qa of 34| can be compared to Q through 
; Qa = z Q 



zok——Ki 
z k ± =K ± 

In the limit v — > there are no incoming or outgoing modes at the horizon. In that case one may require 



that the solution be normalizable [58|. In the discussion following equation (j26|) we show that the v — > 
limit is continuous. 
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momentum scale. Since in this work we are interested in the high momentum asymptotics 
of the solution, we will make sure that the WKB approximation is valid for large Z . As Z 
grows smaller we will inevitably reach a regime where this approximation will break down. 

To apply the WKB approximation to ( fl~4|) we must first bring it to Schrodinger form. 
This will usually be possible for an asymptotically AdS geometry: consider an equation of 
the form 

(V AdS + V')<p = J, (16) 

where T>Ads — d z — 2>Z~ x dz — 1 such that DAds4> = is the Fourier transformed equation 
of motion of a scalar field in AdS^ in appropriate coordinates (see the derivation leading to 
equation (fl4l). The second term on the left hand side of ( fTBT) is given by V = a2(Z)d z + 
ai(Z)dz + ao(Z) and corresponds to a deviation from empty AdS which vanishes at the AdS 
boundary Z = 0. J is a source term. Equation (fl4")) is of this form with 

a 2 {Z) = g{Z) - 1, ai (Z) = g\Z) - ^{g(Z) - 1), (17) 

a (Z) = a 2 (-^ - l) , J(Z) = Z e -*-«zm (lg) 



To bring equation (Tl6l) to Schrodinger form we define 4>{Z) = -j^ip(Z) with b(Z) 



Z 3 exp J Z x ^°fa^ dx. For small values of Z the integral evaluates to —3 In Z which allows 
us to set 6(0) = 1. Next we divide equation ( Tl6l) by yj &py(l + a 2)- The equation of motion 
obtained in this way is a Schrodinger equation, ip" + Vip = 0, with a potential 

v(z) _ zm + (-1 + a ° z ^ 9fl2 + 6Z ( Qi ~ a 'J + z2a i( 2a 2 - go m 

1 ' ~ l + a 2 4Z 2 (l + a 2 ) 2 ' 

One can check that as long as a' 2 (0) = a[(0) = 0, the Schrodinger potential approaches 
VAds(Z) = — 1 — -r§2 in the Z — > limit. The corresponding Schrodinger equation may be 
solved explicitly in terms of Bessel functions (see (I24"I) below). 
For the case at hand we find 

:i + « 2 )(T) 4 41V1 6 



i_5_ _ ^_2zzZ_ ' w w 

422 " ' 



^) = -T7T, ~ — + , v 7 , n2 (2°) 



and 6(Z) = 1 - (j^j . One can see explicitly that as one takes the BH temperature to zero 
(Z — > oo) or considers small values of Z, the potential asymptotes to Vacis- 
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Schrodinger potential for the dilaton 



z min| z 




ct=0.4 



a=2 



FIG. 3: A plot of the Schrodinger potential (equation (|20p ) corresponding to the homogeneous 
term in the equation of motion for the dilaton (equation (|14|) ) for Zq = 40. As a increases from 
to 1, the minimum of the potential increases such that both extrema of the potential vanish at 
a = 1. For a > 1, the zero of the potential gets closer to the origin, decreasing the parametric 
length of the flat region. 



To get a handle on the functional form of the solutions to the Schrodinger equation denned 
by the potential V$ we need to study how V$ depends on a, see figure [31 For a < 1 and 
Z G (0, Zq), one finds that it has two extrema which are given by the roots of a polynomial 
of degree 12. These may be approximated for various values of a. 

As the potential exits the asymptotically AdS region (valid at small values of Z) , it reaches 
a maximum at Z max = ^ (j 1 ^) 1 ^ Z% /3 + O(Z 2/3 ) where V${Z max ) = -1 + 0(Zq A/s ) 
after which the potential is approximately constant until it drops sharply to a minimum at 
z min = (j^) z o + O(Z l ) where it reaches a value of V^(Z min ) = -\ {l+ ^ )2 + O(Z 2 ). 
Here we have implicitly assumed that a ^> Zq 1 but (1 — a 2 ) -1 / 6 is small compared to Z^ 3 . 
Other ranges of a may be treated in a similar manner. After passing the minimum, the 
potential asymptotes to +oo, crossing zero at Z — (1 + a 2 ) -1 / 4 Z + 0(Zq 1 ). With a slight 
abuse of notation, we shall define Z = (1 + a 2 )~ 1//4 Zo. 

For a > 1 the potential has no extrema. As a increases from to 1 the minimum of the 
potential, V$(Z min ), grows closer to —1 such that when a = 1 the potential is almost flat at 
intermediate values of Z; V$(Z) = — 1 — ^fZ~ 2 + ( 4Z ^ zS j . As one increases a further the 
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value of Z becomes smaller, such that the region where the potential may be approximated 
as flat becomes shorter. 

Since a plays an important role in the shape of the Schrodinger potential, we would like 
to make an aside to discuss the values of a in momentum space. Specifically, we will be 
studying the value of (C?$) in the k-, k± plane. Let tan# = k±/k_. This implies that 

v 2 

tan 2 6> = v 2 - 1 + — , 

cr 

so that < a 2 < Large values of a correspond to small 6 (i.e., k± < k-). In order for 

our large momentum approximation to be valid we will see that a should be restricted by 

2/3 

a < Z . If we wish to probe the small 6 region of momentum space at highly relativistic 
quark velocities then we will be restricted to large values of k-. 

Since the potential V$ is technically difficult to deal with, we shall analyze it separately 
in three different regions, using a different approximation in each one (this was also done in 



a somewhat different context in 



41] ). In the near boundary region we use 

v 1 (z) = -i- 



4Z 2 ' 

The near horizon region may be approximated by 

_ \Zlo? + 1 Z (\a 2 + 1) + I 
31 } A(Z -Z) 2 A(Z -Z) 

and the intermediate region by 

(l + « 2 )( 



, 4 

Z 



W) = ; # 2 



One may approximate the regime of validity of these three potentials by ensuring that 
the correction to \Vi(Z) — V$(Z)| is small. Explicitly, for V\ we require 



i 



-l + « 2 )( — ) «1, (21) 



V2 is a good approximation in the range 

'15 



t « Z, (22a) 
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and 

Z - Z > i-, (22b) 

(this is a conservative estimate where a is assumed to be small) and for V3 we find 

Z = C?(Z ). (23) 

Our next step is to solve the appropriate Schrodinger equation for the potentials V\, V 2 
and V3. In each region i we will find two linearly independent solutions which we will denote 
ip l a and For region 1 we find 

^ a = Z l lH 2 {Z) and ^l = Z^ 2 K 2 {Z), (24) 

where here and in what follows we define I u and K v to be modified Bessel functions of index 
v of the first and second kind. 

At small Z we have <f> ~ Z^ip 1 . A series expansion for small arguments of the modified 
Bessel functions in equation (1241) . h(x) ~ x 2 and K 2 (x) ~ x~ 2 , then tells us that the solution 
which does not correspond to a deformation of the Lagrangian (in this case, the solution 
which vanishes at the AdS boundary) is given by 

For V3 we have the solutions 



vV^W/^ [l^^W^) ' (25a) 



= -yftyfZ^ZK*^ I ^v^^J (a2 + 4 i Z ° 2 + 12 I • (25b) 



To understand the near horizon behavior of ip 3 given in equations (1251) . we rewrite the 
modified Bessel functions in terms of Hypergeometric functions 

Ux) = tv 1 ^ (?) " o*i + 1, ^ Y (26a) 



T(z/+ 1) V27 u 1 V 4 

= ir(i/) (|) Fi (-1/ + 1, £) + ir(-iz) (f)Vi (» + 1. x) ' (26b) 

which is valid provided that z/ = |Z a 7^ 0. Defining the Kruskal-like coordinates Y = 
ln(Zo — Z), we find that (Zq — Z)^ v = e ±! i z o vk - Y _ j s now clear that corresponds to 
a solution with modes which are ingoing at the horizon and K v to a linear combination of 
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ingoing and outgoing modes. Therefore ip^ satisfies the near horizon boundary conditions 
we are looking for. For a motionless quark (a — > 0), there are no ingoing or outgoing modes 
at the horizon. Instead one imposes that the solution is normalizable. This implies that we 
should use the Iq solution, which gives us a nice v — > limit. 

In region 2 we need to use the WKB approximation to solve the equation of motion. 



Ensuring the validity of the WKB approximation, 

Z>2 

Z<Z -2. 



<C 1, implies 



(27) 
(28) 



We should also take into account the crossing from V2 < to V2 > at Z ~ Z if this 
crossing occurs in region 2. 

As long as V(Z) = at Z > Z — 2 (or, in other words, a^/Z^ < 1) the potential V 2 is 
strictly negative. This gives us the solutions 



/ 



* = "71 



1 



z 4 



1/1 



I- St 



-ZFi 



I 1 1 - z 



* = 71 



1/4 



1 ^4 



(29a) 



(29b) 



where 



F X (Z) 



1 1 i- 5 -^V 

' 1 1 4' 2' '4'\ZJ ' \Z 



(30) 



is the Appell hypergeo metric function (see for example 63j) which may also be written 
as a linear combination of Elliptic integrals of the first and third kind. As expected, for 
< Z < Zq this function is real as long as its first argument, Z/Z, is smaller than unity. 
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On the other hand, if ayZo > 1 we find that the solutions in region 2 are of the form 

2\ 1/4 



1 



~ ^7f X/ \ 2 V^ ^ 1,1+,,^,'' i 




e ~ZFi(Z)+ZF 1 



1/3 



z) 



cos (-iZF 1 (Z)+iZF 1 + f) 



V 2 (Z)<0 
and 

M\>\v 2 \ 



\V 2 '\ < \V 2 \ 



V 2 (Z)>0 
and 

iv 2 'i>iv 2 r 



(31a) 



e +ZF 1 {Z)-ZF 1 



A, 
i\ 1/4 



4(l +a 2 )9 /4 
Zoa 4 



1/3 



Hi) 



sin (-iZFx(Z) + iZF x + f) 



V 2 (Z)<0 

and 
|V 2 '|>|V 2 | 



|v 2 '| < |v 2 



y 2 (^)>o 

and , 

iv 2 'i>iv 2 r 



(31b) 



where we define the Z independent constant F\ through 



Fx = F 1 {Z) 



2r(I) : 



F-i 



1 



I 1-1- 

4' '4' (1 + a 2 ) 



and Ai and .Bj are Airy functions (see for example j^). 

In all three regions we have chosen our normalization such that the Wronskian of ip a and 
ip b is -1. 

Our next step is to glue the solutions in each region to one another to obtain a continuous 
and differentiable function, if). We do this once with boundary conditions which vanish at 
the AdS boundary, , and once with boundary conditions which are ingoing at the horizon, 
ip + . The gluing should be done at some intermediate points Z\ 2 and Z 2 3 corresponding to 
the transition from V\ to V 2 and V 2 to V3 respectively, such that in each region i the solution 
is a linear combination of the form if)± = A l a± if) l a + A\ ± %p l b (no sum on i, a, b). 

We should choose Z\ 2 and Z23 in the region in which there is overlap between the ranges 
of validity of the potentials V\ and V 2 , and V 2 and V3. Of course our final result will not 
depend on the exact value of Z± 2 and Z 23 . We start with a < 1. Considering Vi, we require 
Z12 ^ v/^f an d ( — 1 + a 2 ) 1//4 Zi2 <C Z . The first constraint comes from requiring that there 
be an overlap with region 2 (see equations (I22al) and (1271) ). The second inequality arises by 
requiring that V\ be a good approximation to V$ (equation (|2ip ). Since the maximum of V$ 



16 



is at Z max with Vq>(Z max ) = — 1 + 0(Zq 4 ^ 3 ) and since V\ asymptotes to —1 at large Z, we 
shall define Z X2 = Z^ 3 ~ Z max . For region 2 we use Z 2 % = Z — 2 due to equation (1281) which 
is compatible with (|22bl) since Z is large. As long as a is small we find that region 1 and 
region 2 have good overlap. To quantify this statement, we observe that V\ — — 1 — and 
for small Z, V2 — > — 1 + (a 2 — 1) y~§^j ■ Our approximation will be good as long as the flat 
intermediate region is captured by these asymptotics. When a > 1, the point Z where the 
potential is most flat is given by the solution to V"(Z) = 0, Z = Z max + O ^Zq 2 ^ 3 ^J. Since 
V{Z) ~ V(Z^ 3 ) ~ — 1 + (a 2 — ^)Zq~ 4 ^ 3 , requiring that a < Z^ 3 will ensure that there will 
be an overlap between regions 1 and 2 for a > 1 as well. 

We proceed to evaluate the various A's. Consider for example the requirement that ip- 
be continuous at Zi_n 

Similar equations hold for the first derivatives of ^_ and two more such sets for the ip + 
solution. Together with Af + = A\_ = and A\ + = A\_ = 1 we have eight sets of equations 
for the eight remaining coefficients. 

The actual values of the A's are listed in appendix [A] Some useful properties of these 
coefficients are 5 

(32a) 
(32b) 
(32c) 

3 _ ~ e*»-». (32d) 

The first three expressions, equations (I32cl) . (I32ap and (132bl) . come about as a result of 
gluing exponentially increasing and exponentially decreasing solutions: an exponentially 
decreasing solution must be enhanced by an exponential factor in order to be glued to an 
exponentially increasing one. Regarding A\_ in equation ( I32dl) . since is exponentially 



4 1 




A 2 


~ e~ 2Z *Al + , 


4-' 


~ e~ 2Zl2 A 2 




~ e Z23+ -, A 



decreasing near the AdS boundary (see the complex valued Gamma functions in fl26b[) ) one 



The actual values of A%_ , A l a+ and A^, are not relevant for the calculation we have in mind and it may be 
that they are exactly zero. Here we take a more conservative approach which allows for an exponentially 
damped coefficient. 
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must boost it when gluing it to an exponential solution in region 2. Similarly, we will find 
that A z a _ ~ e Z23 ~-. 

Using (fT5l) and going back to our less explicit notation (equation (TIE]) ), we find that the 
solution to ( [TBI) is given by 



Z /3 



w(z') v b(z 



-MZ')dZ' 



Za J(Z') 



Z' 3 



z W(Z') v b(Z 



-^ + (Z')dZ'. (33) 



with 



W(Z) 



A 1 



b(Z) 



-A 2 A 2 4- A 2 A 2 



(34) 



< z < z l2 
12 12 z 12 <z<z 23 - 

Z23 < Z < Zq 

Neglecting exponentially suppressed terms, one finds after some algebra that for Z < Z 12 , 
equation ( |33l) reduces to 



-A 3 



b(Z) 



-Z X ' 2 K* 



(Z) / - I 2 (x)dx 

./n X 



z^hiz) 



Zl2 



j(x)^~: 



x 



K 2 (x)dx'. (35) 



Expanding the right hand side of equation (1351) in a small Z power series, and recalling that 
we have defined 6(0) = 1, one finds 



^(^(tO-(M + ^ (0 ))z + j(o)(...)z- 

_ z i/2(l Z 2 + ± z A f Zi2 mVM 

V8 96 JJ Z x 



K 2 {x)dx + 0{Z 7/2 ). (36) 



The last integral may be evaluated by noting that K 2 is exponentially suppressed when its 
argument is large 57]. Therefore, as long as J\/b is not singular in the region of integration, 
one may probe the leading behavior of equation (135]) by expanding it in a power series and 
extending the region of integration to infinity. This will give us 

4> = -\j(q)z 2 - ( -Lj(o)i/(o) + \ J'(o) ) z 3 



24 



-CA;- 2 - -1/(0) J'(0) + J(0)(. ..))Z A - - J(0)Z 4 In -Z + C(Z 5 ), (37) 
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where we have used 



lim J x~ l K 2 (x)dx = + ^2 In QxJ + 2 7 - 1 j , (38) 

r°° 2 1 

lim J x°K 2 (x)dx = - - -it, (39) 

lim/ xK 2 (x)dx = -2 In (-) + 1 - 2 7 , (40) 
x^o _/ x V 2 / 

POO 

/ x n K 2 (x)dx = C n , forn>2, (41) 



where the C n 's are real numbers. For example 

C 2 = ^r, C 3 = 8. (42) 
The constant C in equation ( |37l) is given by the Z independent term in 

p r Mvm Mx)dx (43) 

where corrections of order e~ Zl2 are neglected. 

It is interesting to note that, if J(0) 7^ 0, divergent terms in the z — > limit of z~ 3 d z & = 
k 2 Z~ 3 dz<p{Z) are, due to the logarithm in equation (137|) . non contact terms and a more 
careful renormalization is required jo^l . [o^ ]. This may have some impact on the mesonic 
solutions found in 28J] where the string lies entirely along the boundary. 

It should be emphasized that equations ( |35|) and ( 1431) have been obtained via the WKB 
approximation which was shown to be valid as long as Zq » 1. As we discussed earlier, 
when considering progressively smaller corrections to (Oq>) the WKB approximation will 
eventually break down rendering fT43|) invalid. 

For the case at hand we have J(0) = 0, J'{0) = k^ 1 and b'(0) = 0. Going back to our 
original variable z and original function $, we find 

K ' 3 8 

By our assumptions we find that the leading contribution to the integral fT43|) in the kzQ — > 00 
limit is 

2 1 ~ 

K 2 {x)dx = -7ik. (44) 

k z 2 

As discussed earlier, the z independent term, —-^nk, is the contribution of the near field 
value and is not related to dissipative dynamics (see equation fTTBl and the analysis which 
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follows it). Using the definition of C via (143]) in equation (1371) and 6'(0) = 0, allows us to 
write 

1 1 7 ~ k [ {Z °~ k)2/i / x ( J(x)~kJb(x) 

(£>$) = -—A; + — / K 2 (x) — v v 7 - 1 da; (45) 



V / ^M 16 8tt 

where it is understood that terms of order e~^ z ° k ^ 2/i should be consistently neglected. 
Inserting J(x) = S- e ~ lk -^ x /^ and 6=1— C^) , we find 



(£> $ ) = + A / tf 2 Wl- (4, e-'^-lUx. (46) 



v/Ay M (l-^) 16 87ry V V^ofc/ 

The square root in equation (I46p may be expanded in a Taylor series with an appropriate 
remainder which will give an indication of the validity of the expansion. The oscillatory term 
in the integral may also be expanded in a power series as long as its period is large enough, 
so that the exponential behavior of K 2 enforces the product K 2 e~ tk ~^ to be damped. These 
requirements imply 

• 2 3 < - and Z, > 10. (47) 



3Z% 



Since we are working in the a < Z^ 3 and large Z regime, this expansion is always valid. 
Using equation one may approximate the integral in equation f|46|) to obtain equation 
(CQ) as advertised. 
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APPENDIX A: CONTINUITY AND DIFFERENTIABILITY OF THE 

SOLUTIONS 

In section [Til we have discussed how the gluing of the solutions i/j l a * b , valid in various ranges 
of the Z coordinate, may be parameterized by the sets {^4_} and {^4+}. In what follows we 
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give explicitly the large Z asymptotics of {A±} when a < Z^ 3 . We separate our analysis 
into two regions. Consider first values of a which are small relative to y/Z~o, ay/Zo = e < 1, 
we find 

A\_ = JLz^ s e^ + (Al) 
At = \ (A2) 



3 _ 25 z Q K(-l)-2V2^ , 11 2Z 2/3 -Z K(-l) , n( 2\ /AoN 

~ 32Z e 16Z 4/3 1 J 1 J 

97c; , /q , p Z K(-l) 

3 ZfO -■Z Ji:(-l)+2Z^ + 2 v '2Zo , £ 

512vrZ 7/J 7T 

25i =e ZoX(-i)-2^ _ J^_ e 2zr-ZoK(-i) ) e + (e 2 ) ( A I ) 



and 



64v% 32Z 5/6 



A 2 = e Z K(-l) + O(e*) (J__ % ^M e -2Vm + 0(e 2^ (A5) 

oz-\/7rZo 2 

4I _ 25 2v^Zo-Zo/C(-l) _ 11 Z K(-1)-2Z 2/3 

a+ ~ 32Z 6 16 ^/3 6 

+ ( -^e-^ 4 ' + 275 T 1 /^ y( " l)Z °' 2Z ° 2/3 " 2 ^ e + 0(e 2 ) (A7) 
V 2 1024Z 11/6 y 

x _ -K(-l)Z +2Zl /i +2^m V 

512ttZ 7/3 

/ 25? Zo K{-l)-2Sm 2Z 2 /3 -Z K(-l) \ , / 2X / Ag N 

Here K(—l) is the complete elliptic integral of the first kind [5^], K(—l) = V^ r(3/4) ~ ^ ^ 

— 1/2 

One may check that the aZ = e — > limit coincides with a motionless quark as it should 
(recall that a oc v ). In fact, the analysis was somewhat simplified since the large Z behavior 
of the solutions, ip\ at the gluing points Z12 and Z23 has been expanded in a power series 
in e around the v = 0, large Z behavior. 

— 1/2 

In contrast, the large a case, a > Z , is somewhat more difficult since a large Z 
expansion requires us to know the asymptotic behavior of complex index Bessel functions 
where both the absolute value of the index and the argument is large. Considering equations 
f l26|) this reduces to finding the asymptotic expansion of the hypergeometric function Fi 
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for large values of both its arguments. Finding the asymptotic behavior of hypergeometric 
functions is a notoriously difficult problem. Luckily, the analysis is simplified in the limits 
we are considering and the relevant asymptotic expansion can be carried out, see appendix 

~~ 1 /3 

iBl As long as a > Z we find 

K- = i ^e-»"-, At = 4, (A9) 



7T 



,ZFi- 



TV Zqoc 



A 3 a _ = - - " , At=Oe^ + ^, (A10) 



7T 

A 2 a+ = Oe^^\ < = e-^ + ^, (All) 

^ + = A i + = 0e **+=3*. (A12) 

16Z 7 

where O is a complex oscillating term of order one, which differs slightly in each of the 
expressions. 

To use the approximation in appendix [B] (equation (IB If) ) we must restrict ourselves to 

^ In A. 

a ^> Z . Hence, we are missing the coefficients {A±} for Z 2 ~ a. Nonetheless, since 
the potential is smooth in a and exhibits no distinctive features for this region of parameter 
space, then from continuity and the general arguments given in section [IT] we do not expect 
the A's to behave differently from (|32|) . 



APPENDIX B: AN APPROXIMATION FOR Fi 



Consider the hypergeometric function Fi(b,x) defined by 
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with (b) n = b(b + 1) ... (b + n — 1). We wish to find an approximation for Fi((3N, xN) when 
iV is large. Expanding in a power series and resuming, we find 



F 1 ({3N,xN) = J2 

n=0 

oo 

= E 



1 z n iV n 

1 2 n iV r ' 



^ 7l! ^JV" _|_ 



+ £>(A^ 2 ) 



E 

n=0 



1 x n iV n 



n! /3 n N r ' 



1 - 



1 n(n — 1) 

N(3 2 



+ £>(A^ 2 ) 



1 - 



1 1 fx 



O (0N) 



-2 



Thus 

is a good approximation as long as 
parameters 

ia 



(Bl) 
(B2) 
(B3) 
(B4) 



N(3 2 \/3 



< 1. In appendix |X] we apply this expansion with 



N(3 3 



iV = Z , x~(Z -Z)(« 2 + 4), 



(B5) 



1/3 

(see equations (125]) and (|26l) ). For Z — Z ~ 1, we find that we must require both a>2 
and a Z in order for the approximation to be valid. 
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